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We present calculations of the charming-penguin long-distance contributions to B — * Ktt decays 
due to intermediate charmed meson states. Our calculation is based on the Chiral Effective La- 
grangean for light and heavy mesons, corrected for the hard pion and kaon momenta. We find that 
the charming-penguin contributions increase significantly the B — > Kn decay rates in comparison 
with the short-distance contributions, giving results in better agreement with experimental data. 



I. INTRODUCTION 



Two-body charmless B decays are studied as a mean to detect direct CP-violation and to determine the CKM mixing 
parameters in the Standard Model jl] . In these processes the short-distance contributions can be computed by using 
the factorization approximation for the local operators in the effective non-leptonic Hamiltonian 7Y c ff as argued by 
Bjorken on the basis of color-transparency [|| and recently justified in j|. However non-leading 0(l/rrib) effects 
which appears in the penguin matrix elements may be numerically important. In some cases where experimental data 
are available this approach does meet difficulties (for a discussion see e.g. [Q). An example is offered by the decay 
B — » Kn. The branching ratios for these processes have been measured by the CLEO collaboration ||: 

B(B+ -> K Q n+) = (18.2±to ± 1-6) x 1(T 6 , 

B(B° -> K+TT-) = (17.2±1;| ± 1.2) x lCr 6 , (1) 



and, more recently, by the BaBar B and Belle M collaborations. BaBar gives a value of B(B° — > K + tt ) = 
(12.5t^+i ; 3) x iQ-6. Belle gwes^B{B° -» K+tt-) = (17.4j£J ± 3.4) x 1CT 6 and B{B+ -> K°tt+) = (16.6±?;|±|;J) x 
10 -6 . If one now evaluates the amplitudes for these decays by including only tree and penguin operators, without 
taking into account the charm quark loop, the factorization approximation produces branching ratios too small as 
compared to the data Q . 

To get a better agreement one is forced to include the so-called charming penguin operators, i.e. those operators 
which, being proportional to the large CKM factor V* b V cs and the large Wilson coefficient o% w 1, are not suppressed 
like the tree and penguin contributions. Two approaches can be followed to include these effects. One may define 
effective Wilson coefficients by considering the effect of a charm quark loop treated perturbatively [^o|-p^|; in this 
way also an absorptive part of the non-leptonic decay amplitude is generated. This approach, used together with the 
factorization approximation, seems to produce decay rates in agreement with data, at least qualitatively, as shown 
in some previous papers |l4| [TtJ . In these works the inclusion of the charm quark loop increases the effective Wilson 
coefficients of the strong penguin operators by about 30%, thereby producing B — > Kir decay rates closer to the data. 
More recently char m qua rk effects computed by this method have been included in works dealing with the validity of 
the factorization iggnj . 

A different approach can however be followed. It assumes that the charm quark contributions are basically long- 
distance effects that can be taken into account by including rescattering processes such as, e.g. B — > DD S — > Kw. 
These contributions, first discussed to our knowledge in J2p[| , have been more recently stressed by [Q, where they are 
called charming penguin terms, a nomenclature we shall adopt here. 

The aim of this paper is to present an evaluation of the charming penguins in the B — ► Kit decays going beyond 
the parametrizations of Q. As a matter of fact, instead of considering only D,D S intermediate states, we will also 
consider charmed vector mesons ^| Furthermore, by using phenomenological information from semileptonic decays and 
Chiral Effective Lagrangean, we shall estimate both the real and the imaginary part of the charming penguins. This 
result, being an improved determination of the strong phase, might be of some utility in connections with strategies 
to determine the angle 7 of the unitarity triangle from future more precise data. 

The long-distance absorptive part is essentially due to the rescattering effects of the processes 

£>(*)£)(*) _» Kn. The 

situation is similar to the B s — > 77 decay for which the absorptive part obtained in |22| is comparable to the short- 
distance contribution. As noted above, the dispersive part of the charming penguins, in the present approximation, 
was computed previously in |20| for a number of charmless B decays to two pseudo-scalar mesons as well as one 
pseudo-scalar and one vector meson in the final states. In this paper we present a new calculation of the long-distance 
contribution using the same Cottingham formula p3|j , and more recent information on the semileptonic decay form 
factors. 



1 After having completed this work new data from the BaBar || and Belle |j| Collaborations have been published: BaBar gives 
B(B+ -> K°ty + ) = (18.2+f ^+2 ^) x 1(T 6 and B(B° -» K + n~) = (16.7 ± 1.6±\*) x 10~ 6 . Belle gives B(B+ -> K°ir+) = 
(1-311q.46 ± 0.26) x 1CT 5 and B(B° -> K+ir') = (1.87±g;|f, ± 0.16) x 10~ 5 . Within the errors they are compatible with 
previously released data. 

2 The transition B — > X c is saturated by X c — D + D* at zero recoil |2lJ. Though we are not in this limit, the approximation 
can be used as a guideline. 
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The plan of the paper is as follows. In the next section we discuss the separation between the short-distance 
and the long-distance contributions. In section 3, the absorptive part is shown in terms of the B — ► D^D^ 
and D^Dg — * Kit amplitudes obtained by the short-distance non-leptonic Hamiltonian and the B — > and 
£)(*) — ► A"7r scmilcptonic decay form factors. Section 4 is devoted to the calculation of the dispersive part (the real 
part) of the B — > Kir amplitudes. Finally in section 5 we compare our results to the data and draw our conclusions. 



II. SHORT-DISTANCE AND LONG-DISTANCE NON-LEPTONIC WEAK MATRIX ELEMENTS 



The non-leptonic B — > Kit decay amplitude is obtained by considering the matrix element 

A Klx = < K(p K )7T(p 7r )\in c a\B(p B ) > ■ 



(2) 



The effective Hamiltonian for non-leptonic B decays is the sum of 4-quark tree-level and penguin operators and is 
given by 



T~CcS = 



7f 



V: b V us (ciO¥ + c 2 O u 2 ) + VaVcfaOl + c 2 O c 2 ) - V* b V u c *°* + c a°^j 



(3) 



where Cj are the Wilson coefficients evaluated at the normalization scale /i = m b |ri|,^4,^5| and next-to-leading QCD 
radiative corrections are included. 0\ and 2 are the usual tree-level operators, Oi[i = 3, ...,10) are the penguin 
operators and O g is the chromomagnetic gluon operator. The Cj in eq. (||) are as follows [M : c 2 = 1-105, c\ = 
-0.228, c 3 = 0.013, c 4 = -0.029, c 5 = 0.009, c 6 = -0.033, c 7 /a = 0.005, c 8 /a = 0.060, c 9 /a = -1.283, c w /a = 
0.266. 

In the calculation of the B — > Kit decay amplitude Ak-k we can separate the short-distance and the long-distance 
contributions: 



Artt = A, 



SD 



A 



LD 



(4) 



The short-distance part of the amplitude A$d arises from the operators in (|3|) that give non vanishing contributions 
in the factorization approximation, i.e. Of (i = 1, 2) and Oi (i — 3, 10). In this approximation it is given by: 



A SD (B+ -> K\+) = %f K F B ^(m 2 K ) (ml - ml) x 



V2 



0-4 — 2 ai0 



2 m\ 



(m b - m d ) (m d + m s ) 



A SD (B° -> K+n-) =-^fK F B ^(m%) (m| - ml) x 



V* b V us a 2 - V t * b V ts a 4 + a w + 



2 m\ (a 6 + a 8 ) 



(m fe - m u ) (m u + m s ) 



(5) 



(6) 



where dj = Cj + — — (i=odd) and aj = Cj H — - — - (i=even). 



Numerically, for \V ub \ = 0.0038, V us 
F B ^{m 2 K ) = 0.37, we get 



0.22, Va - 1, H 5 



0.040 and 7 = - arg (V ub ) = 54.8 



^sd(S+ -> i^°7r + ) = 2.43 x 10 -8 GeV 
-4sd(B° -> A+tt") = (1.86 - 1 0.95) x 10~ 8 GeV 



and 



(7) 



As discussed in the introduction by this contribution alone we would obtain branching ratios roughly one order of 
magnitude smaller than the experimental findings. Therefore one has to relax some of the hypotheses. The possibility 
we explore here is to consider non-factorizable contributions. In this case also the operators Of (i = 1,2) are effective; 
these terms will be treated as long-distance contributions, i.e. we will go beyond vacuum saturation and consider 
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intermediate low energy hadronic states in the product of the weak currents in the operators of (^). Also other 
operators in have long-distance contributions, but clearly A^£> is dominated by 0\ due to the enhancement of the 
CKM factor V* b V cs and the Wilson coefficient 02. Therefore we can write: 



A LD = A LD (B+ -> K°7T+) = A LD (B°^K+7r-) = 
Gf „♦„ f d 4 q 



-6(q 2 + n 2 )T(q,p B ,PK,Pir) 



(8) 



where [i is a cut-off separating long-distance and short-distance contributions, while the amplitude T(q,pB,PK,PTr) 
<7 Ml/ Tfj, v is obtained by: 



(9) 



with J M = 67^(1 — 75)0 and J„ — cy v (l — 75)5. Ald has the same value for the two channels B + — > K°ir + , B° 
K + tt~ . We will saturate the T— product of the two currents by inserting D, D* intermediate states, see fig. [lj. 




D,D* 




FIG. 1. The diagram corresponding to the hadronic tensor T M „ in equation (tm. The boxes represent weak couplings. 



To compute the diagram of fig. [l] one has to model the B — > £)W and the — ► Kir weak transitions; this will 
be done following the Chiral Effective Lagrangean approach of [^7j . We list here the effective vertices and currents of 
this theory. 



1) Strong coupling among two heavy and one light meson: 
where 



1+ i> 



{P^r - Pal,) , H a = 7°^7° 



(10) 



(11) 



and 



.M 

1 ''''' ^ l ~r 



(12) 



Here v is the heavy meson velocity, P a ,P* M are the annihilation operators of heavy pseudo-scalar and vector 
mesons made up by a heavy quark and a light antiquark of flavour a (a = 1,2,3 for u, d, s); M. is the usual 
3x3 matrix comprising the octet of pseudo-Goldstone bosons; / is the pseudo-Goldstone bosons decay constant 
130 MeV). 
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2) Strong coupling among two heavy and two light mesons: 



£ Mm = ^TrH a H b v» {£}d^ + id^) ba . (13) 

3) Weak coupling of a heavy meson to pseudo-Goldstone bosons by a g Q 7 M (l — lb)Q current: 

L^ a = %Try„( l-75)#4a- ( 14 ) 



a is related to the heavy meson leptonic decay constant by the formula a = fu Jthh-, valid in the infinite quark 
mass limit. Eq.jl^) generates, for example, weak couplings of D, D* to hadronic final states with 0,1,2,... 
pseudo-Goldstone bosons. 

4) The weak matrix elements < (D, D*)\ J^\B > are parameterized as in [ p7[ , i.e. in the infinite heavy quark limit 
and introducing the Isgur-Wise function for which we use the simple expression 



€iw(v-v') = ^w(w) = (yT^) ' 



(15) 



where v, v' are heavy meson velocities. 



5) The weak matrix elements < K%\ J^\D > and < Kir\ J^\D* > can be computed by the rules given above, i.e. 
using the model of p7[ . It amounts to consider polar diagrams as well as direct production of light mesons in 
the framework of the Chiral Effective Lagrangean This corresponds to the evaluation of the diagrams in fig. 



3 We do not include graphs where the pion is emitted at the first vertex; either they are already taken into account by the 
contributions discussed so far, and one must omit them to avoid double counting, or they are negligibly small, according to the 
estimates in IpOf. 
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V 




(e) (f) 



FIG. 2. The relevant Feynman diagrams used to compute T M „. The box and the circle refer, respectively, to the weak and 
to the strong interaction. Fig. (b) corresponds to three diagrams: the first one has the D state on the horizontal internal line 
and a D* state on the vertical line; the second one presents a D* state on both lines; the third graph has, respectively, a D* 
and a D state. The same rule applies to the other figures. 



Strictly speaking, the evaluation of these diagrams by the Chiral Effective Lagrangean is valid only for soft light 
pseudo-scalar mesons. Therefore, in order to use it in the present context, we have to account for the high momenta 
of the outgoing light mesons. To this aim we introduce two modifications: i) we keep the full propagator in the pole 
contributions instead of the expressions in the soft pion limit (a similar use of the full D* propagator to go beyond 
the soft pion result has also been given in p£|) ; ii) we introduce a form factor in the strong coupling constant of 
light and heavy mesons (a similar approach is used in semileptonic decays [p9[). Let us consider the D*Dir coupling 
that can be written in general as follows: 

(D(p / )7r(jv)|I>*( P ,e)) = G D *nn(e-p„). (16) 



From @ we have 



_ 2m D g 

<^D*Dtt — 7 j U'j 

Jit 
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where g is independent of the heavy flavor and is predicted to have the value (see p{| and references therein): 

g re 0.40 . (18) 

We shall neglect I/tuq effects in this case and shall adopt this value. There is, however, an important point to be 
discussed here. The expression ([l6]) is derived from the low momentum Chiral Effective Lagrangean ( |To| ) . The neglect 
of higher pseudo-Goldstone boson derivatives is however not justified in our case, as the pion and kaon have momenta 
of the order of ms/2. To include this effect we modify ( |ri| ) as follows: 



GW = ^-F(\p n \) , (19) 

where F(|p w |) is a form factor normalized as F(0) = 1. This form factor can be evaluated by using the constituent 
quark model. For this purpose, let us introduce the heavy meson wavefunction in the momentum space: ipi)(k), 
where k = = |<fi — gal/2 is one half of the relative momentum of the two component quarks (whose momenta are 
respectively qi and 9*2)- The coupling constant g in the soft pion limit is proportional to the overlap of the D* and D 
wave functions: 

g cx J d 3 kr D (\k\)^D*(\k\)f(k) (20) 

where f{k) is some smooth function whose precise shape depends on the particular model one employs^, but it is 
irrelevant for our purposes. We therefore get, if the pion momentum is p v 



nm = 



d 3 k ij)* L 



Ptt_ 

2 



d 3 kr D (\k\)^(\k\)f(k) 



^D*(\k\)f(k) 

(21) 



where the denominator has been introduced to normalize correctly the form factor. We employ for the wave function 
the expression^: 

^(fe) = ^.(fe)«e-* (22) 



which corresponds to an average quark momentum 

< k >= - (23) 
a 

inside the meson and we assume a constant value for the smooth function f(k). We plot in fig. || the value of the form 
factor as a function of the pion three-momentum for two values of a, corresponding, respectively to < k >= 400 MeV 
(upper curve) and < k >— 300 MeV (lower curve). 



4 For a particular calculation using this approach see [ pp[ . 

5 This expression is sufficiently general; it also corresponds to fits of the wave function in particular constituent quark models, 
see e.g. HS. 
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FIG. 3. Form factor F(p) for the D*Dtt strong coupling as a function of the pion momentum p — in GeV. Normalization 
is for the soft pion limit. The two curves correspond to < k >= 400 MeV (upper) and < k >= 300 MeV (lower). 



For \pn I ~ itib /2 we get therefore 

= 0.065 ± 0.035 , (24) 

where the central value corresponds to < k >— 350 MeV and the higher (lower) to < k >— 400 (< k >= 300) MeV. 

It should be observed that also the Lagrangean Crhtt-k and the weak current have corrections from terms 
containing higher order derivatives. However these terms do not contribute to the imaginary part of Ald; on the 
other hand hard pion and kaon effects in the calculation of to the real part can be taken into account by an appropriate 
cut-off of the Cottingham formula. We will discuss the problem in section 4. 



III. IMAGINARY PART OF THE LONG-DISTANCE CONTRIBUTION 



It can be easily seen that only diagrams ||c and |^d contribute to the discontinuity of Ald- The diagrams ||a and ^> 
have no discontinuity, whereas the diagram ^e vanishes in the chiral limit and we neglect it altogether. As to the 
diagram |]f, both its imaginary and its real part vanish, as it can be easily checked. 

We use the following kinematics: 

p^ = tob« m = (ra.B,0) , = m D v ,fl , q = p-p D M . (25) 

The discontinuity of the diagrams of figs. (|b,|=) gives 

J^[ S +(l 2 - m D s ) S +(PDM ~ m L>) X 

x A(B -» DWDM) A^WflW -> Kit) = 



m D 



16ir 2 rriB 



Vuj* 2 - 1 / dnA(B -> DWDW)A(I)W])W -> Kir) , (26) 



m B + m 2 D — m 2 D ... . . _ — 

where to* = and the angular integration is over the directions of the vector v' = nyuj* — 1. Our 

results are as follows: 



A(B -» D S D) = -K{m B - mu)(l + u*) , 
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A(B -> D*(r), q)D*(e, «')) = Km Ds rj^e* a (ie^t/ V - g» a (l + lu*) + v a i/ ) , (27) 



G i 



where K — —=V* b V cs a2\ZmBm,ofDAiw( 1 ^'*) ■ On the other hand we have 
V2 

A(D s D^K,)=i^ m f}^\ 

JnJK 

-Ptv -Pk +P-k ■ {q-PK){PK ■ q " m 2 K )/m 2 D , 
(m D v' - Ptv) 2 - mjj, 

Arr>*< \rw ,s i v \ ■ @9 F (\P*\)) 2 y/m D1 m D . 

A i D s (V, q)D (e, v') -> Xtt) = i — v - eyr) a x 

JirjK 

ma* vIvk G aX (p w ,p K ,v') 



where 



(m D *v' -Pn) 2 - m 2 D (m D *v' - p^) 2 - m 2 D , J ' 



G aX {p^ PK ,v')= ~(v' ■q)[g a \pK-p,)-plp X K) - (q-Pir)(v' a p X K -g° X (v'-pK) 



(28) 



-q^p^PK-V^-V^pK-p^j . (29) 

We have not written down other amplitudes with no contribution to the discontinuity. Our numerical results obtained 
for the central value in Eq. ( pi|) are reported in table 1. 

TABLE I. Numerical values for the imaginary part of Ald- Units are GeV. First column refers to the D,D S intermediate 
state, the second column to the D* , D* intermediate state. 



D, D s 


d*, d: 


Total 


1.45 x 1(T S 


0.89 x 1CT 8 


2.34 x 1(T S 
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IV. REAL PART OF THE LONG-DISTANCE CONTRIBUTION 



The real part of the diagrams and |^d could be computed by a dispersion relation, from their imaginary parts; 
however this procedure suffers from the uncertainty related to possible subtractions. A way to include them is to 
follow a Feynman diagram approach, using the Effective Lagrangean discussed in section 2. This basically amounts 
to including the real parts of figs. ||a and ||b and also the other diagrams of fig. || which, as we have seen, do not 
contribute to the imaginary part. 

To compute Re Ald we first observe that we can change the integration variable in (||) from q = ps —£>£>(*) to the 
momentum £ defined by the formula 

q = p B - P D M = (m B - m DM )v - £ . (30) 

We note that, by this definition, I measures the virtuality of the intermediate state while the velocities of the two 
hadrons coincide; one can always make this choice using the reparametrization invariance of the Heavy Quark Effective 
Theory |29[| . The cut-off \xt corresponding to the momentum £ can be evaluated, within our model, as follows. For an 
on-shell meson with momentum pu = ttid v, the two constituent quarks have momenta p c = k, p q = —k in the meson 
rest frame. Adding f shifts p c — k — > k + 1; therefore the argument of the wave function appearing in the calculation 
of the Isgur-Wise function, instead of \k\, is \k + £/2\. This corresponds to introducing in the amplitude the form 
factor with F given in fig. ^[ Moreover, we have to implement the condition that the residual momentum 

k + I of the heavy quark does not exceed the chiral symmetry breaking scale, i.e. a mass scale around 1 GeV. Since 
we assume < k >= 350 MeV, this gives the condition \£\ < 0.65 GeV. The smooth form factor F{\£\) (\£\ < 0.65 GeV) 

<.0.65 GoV ^ t 

can be substituted by a sharp form factor 0(fig — \£\), and fj,g can be fixed by imposing / F(x)dx = / dx. 

Jo Jo 

This procedure gives 

fi e ps 0.6 GeV . (31) 

On the other hand £q is of the order £ 2 / (2m c ) and is therefore negligible in the large heavy quark mass limit. Therefore 
we conclude that a cut-off fj,£ on \£\ as given in ( |3~i]) reflects in a cut-off fij on £ 2 . In passing, we observe that higher 
derivative corrections to the Lagrangean Chh-h + Chhiht and the weak current L M produce negligible effects due to 
this cut-off procedure. 

Having fixed the cut-off we now compute the £ integration by performing a Wick rotation: £° — > i£° and changing 
integration variable from \£\ to L 2 = — ((if 1 ) 2 — f 2 ). We get therefore a Cottingham formula [^3| as follows: 

Here f D(t) =< £> W |&7^(1 - 7 5 )c|B >, =< Ktt\ct/^(1 - 7 5 )s|£> ( * ) >■ Both and can be found in §7|f| 

The various contributions to the , corresponding to the different graphs in fig. H, are as follows: 

ftR = - f -fPfe» , (33) 

^ = 1 UfK (pu^-p,) 2 -^ + 

2 gfD*m D , 1 a 1 

2 e M«/37 £ Prr Pd* ' 



h 



U fi< (pd* - Ptt) 2 - m 2 D , 

(g,) _ , 4 g 2 fp*m D , (e ■ p„) 

fw !k [(pd* - Pk) 2 - m%} [(p D * - Ptt- Pk) 2 - m 2 D , 



* M 



As we have already said, we correct the heavy meson propagator to include the hard pion momenta. 
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, \ (PD* -Pk) - {Pit -Pk) , s 
{PK)^ 2 (PD* ~P« ~PK)^ 



D * 



4 g 2 fD*m D * 



U fK [{PD* ~ Pit) 2 ~ m 2 D ,] [{p D . - p„ 



Pic) 



(Po -P^ie-PK) + {Ptv -Pk){s -Ptt) 



(P£>« 'Ptt -p 2 D ,){s-p K ) ~ (PD* -Pk)(£-Piv) 



Pr>* -P-k){p-k - Pk) ~ (pd* ■ Pk)(pd* ■ Pn) 



wo 



.f D *rn D * 



PD- ■ (Pir - Pk) 



* 



/tt /k (PD* 

4 g 2 fp*m D . _ 

U Ik 
(pd* - P-K 



Pit - PkF 



(Ptt +Pk) 



{p-Pv- Pk)^ 



a IS 

taP-iS P" P D * Pk e 



[(j>D* - 
■VK)ii , 



p*Y 



n 2 D ,] [(p 



D* 



PK) 



(34) 



where g = gF(\p„\). The results for fi£ = 0.5 4- 0.7 GeV are reported in table O. 

All the terms in the previous equations containing the factor g are corrected by the form factor f (|j7,r|); the terms 
/j( 2a ) and h^ 2 ^ should also contain their own multiplicative form factors F a and F e (this holds both for D and D* 
intermediate states). We have not written down them explicitly for the following reasons. F e would multiply a term 
which vanishes in the chiral limit and does not affect the final result. On the other hand F a (q 2 ) would represent a 
correction for hard pion and kaon momenta. In the region of high q 2 , i.e. q 2 S [m 2 D , (tub — mn) 2 ], F a (q 2 ) is smooth 
and can be put equal to 1, i.e. to the value corresponding to the soft pion and kaon limit. This result can be inferred 
from the analogous behaviour of the form factor F^^' K (q 2 ) describing the coupling of a heavy and one light meson. 
In the high q 2 region the contribution of the low lying pole B(Q + ) to this form factor vanishes in the chiral limit 
and the form factor is dominated by the direct coupling displayed in Fig. 2a, giving, as a result, the simple formula 
/b//tt- Several numerical analyses confirm a smooth behaviour of F ( f^ 7T (q 2 ). For example in the quark model of |3^] 
v (q 2 ) increases by 40% in the range 15 — 26 GeV 2 ; this behaviour can be fitted by a formula 



^0 



^0 



(q 2 ) oc 1 - 



(m B + to a ) 2 ) 



(35) 



with 771a — 2.5 GeV. In other models a smoother or similar behaviour is found, see the discussion in 
For D decay this formula would hold with m B — > mo ■ 



(m D + toa) 2 ) 



(36) 



However, it would be a too strong assumption to assume that F a (q 2 ) is given by this formula; therefore we put F a = 1 
in the sequel and we shall enlarge the theoretical uncertainty by an extra amount that we estimate ±50%. 

Let us finally observe that, while the numerical value of F a does not vary significantly in the range q 2 € [to 2 ^, {m B — 
Tod) 2 ], its value for q 2 = (777,73 — rnu) 2 is formally suppressed by one power of (I/mi) a s compared to F a {m 2 D ), 
assuming a qualitative behaviour analogous to Fq > ~** (<? 2 )- 



TABLE II. Numerical values for the real part of Ald for fie — 0.5 -j- 0.7 GeV. First column refers to the D intermediate 
state, the second column to the D* one. Units are GeV. 





D 


D* 


Total 


0.5 


-4.66 x 10 _u 


1.62 x 10~ B 


1.15 x 10" s 


0.6 


-7.77 x 10~ u 


2.79 x 10~ s 


2.01 x 10 -8 


0.7 


-1.19 x 10~ 8 


4.40 x 10~ s 


3.21 x 10~ 8 
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V. CONCLUSIONS 



Our numerical results show that the long-distance charming penguin contributions to the decays B — > Kn are 
significant. These results agree qualitatively with a phenomenological analysis of these contributions given in [Q. In 
particular, we found that the absorptive part due to the D, D s states is somewhat bigger than that from the D* , D* 
states, but of the same sign. The real part due to the D*, D* states is however 3 — 4 times bigger and opposite in sign 
to the contributions from the D, D s states. As shown in table 1 and table 2, the total contribution for the real part 
and absorptive part are of the same order of magnitude, at the 10 -8 GeV level. The results for the branching ratios 
are as follows: 

B(B+ K\ + ) = (2.4+H) x 1(T 5 

B(B° -» K+n-) = (1.5±i;f) x 1(T 5 . (37) 



The central values in this equation correspond to the central values of the parameters g, < k > (i.e. i* 1 ^]) and m. The 
theoretical uncertainties in the branching ratios are obtained by varying the parameters in the ranges g = 0.40 ± 0.08, 
< k >= (350 ± 50) GeV and (it = (0.6 ± 0.1) GeV. An extra theoretical uncertainty of 50% has been added in 
quadrature to ReA^D. The results in ( j37[) agree with the experimental values in eq. (0). 

The existence of two different contributions with different strong and weak phases produces a direct CP violation 
in the decay B — > X ± 7r =F . As a matter of fact, we find for 

= T(B ^K-n+)-T(B^K+7r-) 

r(B° -» K-n+) +T(B° -» K+n-) [ ' 



A 



the value A — +0.21 (for 7 = 54.8°), which is compatible with the recent results from CLEO j34j. On the other hand 
the present model produces no CP asymmetries for charged B decays into Kir. 

Let us finally discuss the problem of the scaling of our results with mj,. Assuming, as in J3|, that F^'~ >7 '(0) scales 

Q CD ] Asd scales as ml . On the other hand the charming penguin contribution Ald is non leading due 
m b J 

to the various form factors which vanish for mj,, m c — > 00. Therefore the following conclusion can be drawn. The 
charming penguin contributions violate naive factorization and are of leading order in a s , as they arise from non 
perturbative calculations; nevertheless they do not contradict the results of |^,|35| since they are suppressed in the 
infinite heavy quark mass limits. In spite of that, as mb,m c in reality are finite, the charming penguin numerical 
contribution to the branching ratios in B — > Kir decays is significant, basically due to the enhancement of the CKM 
matrix elements. 

In conclusion, we believe that the charmed resonance contributions we found seem to be capable of producing the 
charming-penguin terms suggested in 0] within theoretical errors. This method could then be applied also to other 
two-body nonleptonic B decay channels, such as the B — > tttt and the B — > Krl . We only mention here a result 
for the charming-penguin contribution in B° — > ir + ir~ decay. This contribution can be obtained from the results we 
obtained for B° -> K+ 7r using SU(3) symmetry for the weak current matrix elements < K + ir \(cs)l\D s > and 
< tt~ -k + \{cd) l\D~ >. Thus, in the SU(3) limit(by ignoring the K — tt mass difference) p7j 

A LD (B° - ^+tt-) = {V cd /V cs )(f K /U) x A LD (B° - K+n-) . (39) 



Since the B — > ir + ir~ decay is dominated by the tree-level operators, this charming-penguin contribution behaves 
like the penguin terms and reduces the B — * tt + it~ decay rate by a small amount. The details of this work will be 
given elsewhere together with the results for B — > Krj ' . 
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